Orbital-selective Mott transitions in the anisotropic 
two-band Hubbard model at finite temperatures 



C. Knecht, N. BliimerQ and P. G. J. van Dongcn 
Institute of Physics, Johannes Gutenberg University, 55099 Mainz, Germany 
(Dated: February 2, 2008) 

The anisotropic degenerate two-orbital Hubbard model is studied within dynamical mean-field 
theory at low temperatures. High-precision calculations on the basis of a refined quantum Monte 
Carlo (QMC) method reveal that two distinct orbital-selective Mott transitions occur for a band- 
width ratio of 2 even in the absence of spin-flip contributions to the Hund exchange. The second 
transition - not seen in earlier studies using QMC, iterative perturbation theory, and exact diago- 
nalization - is clearly exposed in a low-frequency analysis of the self-energy and in local spectra. 
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The Mott-Hubbard metal-insulator transition - a non- 
perturbative correlation phenomenon - has been a sub- 
ject of fundamental interest in solid state theory for 
decades^ Recently, this field became even more exciting 
by the discovery^ of a two-step metal-insulator tran- 
sition in the effective 3-band system Ca2-. T Sr T Ru04, 
for which the name orbital-selective Mott-transition 
(OSMT) was coined^ The Ca2_ x Sr. r Ru04 system was 
investigated theoretically in detail by Anisimov et 
alX within the local density approximation (LDA and 
LDA+U) and within dynamical mean-field theory^ 
(DMFT) solved using the non-crossing approximation 
(NCA). The underlying assumption of a correlation 
(rather than lattice-distortion) induced OSMT found 
support in further band structure calculations^! and 
strong-coupling expansions for the localized electrons in 
the orbital-selective Mott phased 

Microscopic studies of the OSMT usually consider the 
2-band Hubbard model H = Hi + H 2 , where 



{ij)mcr itn 

includes hopping between nearest-neighbor sites i,j with 
amplitude t m for orbital m £ {1,2}, intra- and inter- 
orbital Coulomb repulsion parametrized by U and U' , 
respectively, and Ising-type Hund's exchange coupling; 



for spin o G {t> 4}- I n addition, 



contains spin-flip and pair- hopping terms (with 1 = 2, 
| =1 etc.). In cubic lattices, the Hamiltonian is invariant 
under spin rotation, J z = Jj_ = J; furthermore U' = 
U — 2 J. In the following, we refer to Hi + H 2 in this 
spin-isotropic case as the J-model and to the simplified 
Hamiltonian Hi as the J z -model. 

Liebsch2*i2iii questioned the OSMT scenario for 
Ca2- 2; Sr 2 ;Ru04 on the basis of finite-temperature quan- 
tum Monte Carlo (QMC) calculations (within DMFT) 
for the J z -model using J z = U/4, U' = U/2, and semi- 
elliptic "Bethe" densities of states with a bandwidth ra- 



tio W2/W1 — 2. Additional studies using iterative per- 
turbation theory (IPT)Ai seemed^ to confirm his con- 
clusion of a single Mott transition of both bands at the 
same critical [/-value. Meanwhile, Koga et al. found an 
OSMT using exact diagonalization (ED), applied to the 
full J-modeljiii but not for the J 2 -modelii^ Consequently, 
the OSMT scenario was attributed to spin-flip and pair- 
hopping processes. 

Very recently, four preprints appeared^i&ALiS in 
which the OSMT was investigated in detail within the 
DMFT framework. Ref. Il5l applied the Gutzwiller varia- 
tional approach and ED to the J-model at temperature 
T = and confirmed the existence of an OSMT, provided 
that the ratio W2/W1 of the two band widths is suffi- 
ciently small. Interestingly, Ref. also suggested the 
existence of small spectral weight near the Fermi level of 
the narrow-band subsystem in the orbital-selective Mott 
phase. Similar results were obtained by de' Medici et 
al.f& who used slave-spin mean field theory (which is 
closely related to the Gutzwiller method). Arita and 
Heidi! used the projective QMC method to investigate 
the J-model at T = and demonstrated a first OSMT for 
J = U/4 and U = 2.6 (in units of half the narrow-band 
width). Finally, Koga et alii^ used QMC to characterize 
the OSMT for the J-modcl at finite T on the basis of 
spin, charge and orbital susceptibilities as well as spec- 
tral functions; they further showed that additional hy- 
bridization between the bands smears out the OSMT at 
T = 0. Remarkably, the problem originally investigated 
by Liebsch, 9 : 10 : 11 i. e., the J-model with U' = U - 2J Z 
was addressed only in Ref. |Tfl where de' Medici et al. 
found an OSMT at T = within slave-spin mean-field 
theory, in contradiction to Licbsch's and Koga's earlier 
results. Since the slave-spin method is essentially un- 
controlled, the existence of an OSMT in the J z -modcl 
therefore remains unclear. 

The goal of this work is to clarify this issue and to 
establish whether Hi can be regarded as a minimal model 
for the OSMT. In the following, we will sketch our high- 
precision DMFT algorithm, which supplements QMC by 
a high-frequency expansion of the self-energy, before we 
discuss relevant observables and data-analysis techniques 
and present numerical results. 
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Quantum Monte Carlo scheme - Within the DMFT, 
the lattice problem is mapped onto a single-impurity An- 
derson model with general hybridization, which has to 
be determined self-consistently^ DMFT algorithms are 
characterized by the numerical method employed in the 
solution of the impurity problem and by the iterative 
procedure used in order to establish self-consistency In 
this work, the imaginary-time Green function G(t) of the 
multi-band impurity is obtained in discrctized form from 
Hirsch-Fye QMC simulations. Since the lattice Dyson 
equation is most easily formulated in the frequency do- 
main, Fourier transformations (FTs) are needed twice 
per DMFT self-consistency cycle. Since naive FTs would 
violate the analytic properties of Green functions and 
self-energies, all DMFT-QMC codes use either a special 
transformationi^ or interpolate the discrete QMC data 
by cubic splinesi^S Recently, it was realized2i*2£ that 
natural boundary conditions are not well suited for mod- 
eling imaginary-time Green functions and that misfits 
can be reduced by allowing for non-zero second deriva- 
tives of G at the boundaries r — > 0, r — > However, 
this improvement cannot eliminate misfits arising from 
higher order derivatives of G, which may be large at 0, f3 
and are neglected in cubic splines. 

Our DMFT-QMC method further improves on this sit- 
uation by applying the FTs to a difference Green function 
AG = Gqmc — Ghfe, where the model Green function 
Ghfe is determined from a high-frequency expansion in 
such a way, that the second derivatives of AG(r) vanish 
and the higher derivatives are reduced at the boundaries. 
Consequently, natural cubic splines are appropriate for 
AG. The full FT to Matsubara frequencies is then ob- 
tained as G(iui n ) = AG(iuj n ) + Ghfe(*<-<-V 1 ), where the 
last term is computed directly via the Dyson equation 
from the model self-energy Shfe(*w„). The latter is con- 
structed on the basis of the asymptotic expansionSiiS^ 
S 7 M = E^ 7 U M {n p ) + 4 2) M + 0(4,) with 

E^ 2) (w) = -V U a ^Up^({n a np)-(n a )(n )). 

id * — '",P/7 

Here a, (3 and 7 are multi-indices, combining spin and 
orbital degrees of freedom: a = (m, cr); the interaction 
matrix is defined by U ma . mS = U, U ma< ma' = U'-J z S aa >. 
For the one-band model, whereSi H a 2 \id) = [U 2 (n s ){\ - 
(ng))]/u>, the exceptional accuracy of our DMFT-QMC 
method has already been established in comparisons to 
semi-analytic extrapolated perturbation theory for the 
insulating phased Since the algorithm reduces the dis- 
cretization error to the inevitable Trotter error, reliable 
QMC-data can be obtained already from a rather coarse 
discretization; the following results are based on QMC 
for At = 0.4 unless indicated otherwise. The DMFT 
cycle was typically iterated 20 times with 10 6 sweeps af- 
ter convcrgency, using 1000 Matsubara frequencies. Full 
details of our algorithm will be presented elsewhere^ 

Results - In the following, we present QMC results for 
the J 2 -model with J z = C//4, U' = U/2 and semi-elliptic 
densities of states with full bandwidths W\ = 2, W2 = 4, 
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FIG. 1: (Color online) b) Ratio r = Znarrow/Zwide of the 
discrete QMC estimates of the quasiparticle weights of both 
bands (shown in the inset) versus interaction U for vari- 
ous temperatures T; lines for U > 2.1 represent piecewise 
quadratic fits, a) Same data with linear offset. 

for the "narrow" and "wide" band, respectively. 

A traditional criterion for metal-insulator transitions 
is the quasiparticle weight Z = [1 — rfRe E/c?w| w= o] 1 
[for self-energy £(w)] which in the context of QMC 
simulations is estimated in a secant approximation as 
Z « [1 - ImE(i7rT)/(7rT)]- 1 . Evidently both definitions 
recover the limit Z — > 1 in the absence of interactions 
(when S = 0); however, in the insulating regime the dis- 
crete approximation necessarily remains finite, whereas 
the true Z vanishes exactly. As a consequence, metal- 
insulator transitions are expected to appear washed-out 
at finite temperatures. This is seen in the inset of Fig. 

for T = 1/32: Z drops only by about 60% when 
the narrow band becomes insulating at U w 2.0. A 
second transition is not visible on this scale. The ra- 
tio r = ^narrow /^wido shown in the main panel of Fig. 

is more illuminating since we can clearly distinguish 
three regions: in region I (U < 2.0), r is of order unity 
with a sharp decrease near the boundary; in region II 
(2.0 < U < 2.5), r is nearly constant and of order 0.1, 
and in region III (U > 2.5), r increases with nearly con- 
stant slope (until it approaches a finite limit). In order 
to analyze the boundary between II and III, we have 
performed piecewise quadratic fits to the QMC data for 
both regions separately (lines for U > 2.1 in Fig. 
The kinks at the boundary, barely visible on this scale, 
but resolved after subtracting a linear offset in Fig. [TJi, 
indicate a second transition at least for T < 1/32 (for 
which all QMC data fall on the fit curves). 

The intrinsic ambiguities associated with the discrete 
estimate for Z can be overcome in a full low-frequency 
analysis of the self-energy as shown in Fig. [21 Here, 
the data points represent QMC estimates of the product 
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FIG. 2: (Color online) Low- frequency analysis of the self- 
energy for T = 1/32: QMC data (crosses: for At = 0.40, 
circles: for At — 0.32) for the product of frequency ui and 
self-energy E at the Matsubara frequencies is extrapolated 
by cubic polynomials in ui for the narrow/wide band (up- 
per/lower panel). Inset: extrapolated value for lu — » i0 + . 

wE(w) at the Matsubara frequencies iu> n ; these products 
are real-valued since the self-energy is purely imaginary 
on the imaginary axis due to particle-hole symmetry. The 
lines, given by cubic polynomials in lu, are expected to 
extrapolate to at least linearly within a metallic phase 
(where £ is regular) and to a finite value within an insu- 
lating phase. By this criterion, the narrow band (cf. Fig. 
Et) becomes insulating for U > 2.0 while the wide band 
(cf. Fig. Eb) remains metallic up to U = 2.4. Note that 
the QMC results obtained from discretizations At = 0.4 
(crosses) and At = 0.32 (circles) agree on the scale of 
the figure. The extrapolated product co'E,(u))\ lA} — i0 +, a 
measure of the singularity in the self-energy and roughly 
proportional to the expected gap, is shown at better reso- 
lution as a function of U in the inset of Fig. 03. We clearly 
see the presence of two distinct transitions for the narrow 
and wide band at U = U c i ~ 2.0 and U = U C 2 ~ 2.5, 
respectively. Corresponding results for higher and lower 
temperatures (T = 1/25 and T = 1/40, not shown) are 
barely distinguishable from Fig. except for U « U c \. 

The spectral function N{lu) = — ImG(o;)/7r, obtained 
from QMC by analytic continuation using the maximum- 
entropy method (MEM), is depicted in Fig. [31 Evi- 
dently, N(u = 0) is accurately pinned at its noninteract- 
ing value up to U = 1.8 for the wide band (lower panel) at 
T = 1/40 and nearly pinned for the narrow band (upper 
panel). The narrow band develops a gap for U > 2.05, 
where the wide band is still clearly metallic with a near- 
flat DOS at low frequencies. Only for U = 2.6 a gap ap- 
pears also for the wide band; note that the corresponding 
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FIG. 3: (Color online) Spectral function of narrow/wide band 
(upper/lower panel) from QMC + MEM for T = 1/40. 
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FIG. 4: (Color online) QMC estimates for value N(0) of spec- 
tral function at Fermi energy versus interaction. Full/empty 
symbols correspond to the narrow/wide band at At — 0.4; 
crosses to results for At = 0.32 at T = 1/32. 

curve has a striking similarity to the narrow-band spec- 
trum for U — 2.05. Also noteworthy are the sharp band 
edges at low frequencies (where MEM is most reliable) 
and the rapid decay at high frequencies; apparently, the 
high quality of our data has kept artificial broadening to 
a minimum. The full dependence of N(lo = 0) on interac- 
tion and temperature is better resolved in Fig. which 
clearly exposes the orbital-selective Mott transition: at 
U c i a sudden decay to is observed only for the narrow 
band, while the wide-band value is reduced only by some 
20%. Evidently, the second band becomes insulating only 
for U C 2 ~ 2.5. Both transitions become sharper at lower 
temperatures. A discretization dependence is visible for 
T = 1/32 only at U = 2.0. 

Finally, we discuss the intraorbital double occupancy 
D m = {n m -\n m i). As seen in Fig.Eb, -D w ide barely shows 
any features near U c i , while Z? na rrow is reduced by about 
50% from [7=1.9 to [7 = 2.1. We conclude that the wide 
orbital remains itinerant at a point where the narrow or- 
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FIG. 5: (Color online) Intraorbital double occupancy D m = 
(n m |n m |). Lines indicate piecewise cubic fits. Upper panel: 
a linear offset exposes the features of the wide band. 

bital is already localized. A second transition of the wide 
band is clearly seen only after adding a suitable linear 
term in Fig.|3Ji. Here, both transitions appear as regions 
of enhanced temperature dependence and as kinks. Note 
that up to the discontinuity at U c \ for T = 1 /40 and on 
a broader scale the kink is even more pronounced at U C 2 
than at U c ±. Once again we find two phase transitions, 
each of which can be associated with a Mott transition 
of one orbital. Signatures in observables associated with 
the other orbital reflect the fact that phase transitions 
usually leave traces in every observable. 



Discussion and Outlook - We considered a two-band 
Hubbard model with distinct band widths W2 = 2Wi 
and interaction parameters U,U'~U /2, J z = U/4. Using 
high-precision QMC calculations, in which the correct 
high-frequency behavior of the self-energy is carefully 
taken into account, we showed that this J z -model con- 
tains two successive metal-insulator transitions and can 
hence be considered as a minimal model for the OSMTs, 
observed experimentally in Ca2- a; Sr a ;Ru04. These tran- 
sitions are particularly clearly revealed by the low- 
frequency behavior of the self-energy, but also visible 
in the spectral functions, quasiparticle weights, and in- 
traorbital double occupancies. Our high-precision data 
correct earlier QMC result92*i2*ii by Liebsch. They also 
contradict early ED results^ by Koga et al., who may 
have missed the relatively narrow OSMT region due to 
ED's finite energy resolution. On the other hand, our 
results are in qualitative agreement with de' Medici et 
aVs recent slave-spin mean-field phase diagram for the 
J 2 -modelji& which also contains an orbital-selective Mott 
phase, albeit at slightly larger [/-values (at T = 0). 

Our results imply in particular that isotropy of the 
Hund exchange is not a prerequisite for the existence of 
OSMTs, as was suggested by Koga et al~ on the basis 
of their ED study. In future work it will be of interest 
to investigate the effect of other terms in the Hamilto- 
nian, such as hybridization terms of various symmetry; 
first results in this direction have already been reported 
in Ref. 0. It would also be of interest to extend the 
calculations reported here to lower temperatures, to the 
more realistic 3-band case and to magnetic phases. 
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